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Abstract 

Hydrophobic molecules such as oils and certain drugs can be encapsulated between the two 
leaflets of an amphiphilic bilayer in both lipid and polymer systems. We investigate the case where 
the hydrophobic molecules are incompatible with the amphiphile tails and so form droplets. Using 
a coarse-grained mean-field model (self-consistent field theory, or SCFT), we find that droplets of 
a wide range of sizes have the same characteristic lens shape, and explain this result in terms of 
simple capillarity arguments, consistent with the measured variations of surface concentrations of 
amphiphile in the bilayer and in the monolayers that cover the droplet. We study the effect of 
the strength xbo of the repulsion between the hydrophobic liquid and the amphiphile tails on the 
droplet shape, and find a gradual fiattening of the droplet as xbo is reduced. The droplet remains 
at least metastable even at very low values of xbo- This is in contrast to the behavior as the length 
of the hydrophobic molecules is varied. Specifically, if these molecules are at least as long as the 
amphiphile tails, increasing their length further is found to have little effect on the droplet shape, 
while reducing their length below this value quickly causes the droplet to become unstable. 
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I. INTRODUCTION 



Amphiphilic molecules in solution form bilayers for a wide range of molecular architectures 
and experimental conditions These structures may form both from lipids 3| and from 
block copolymers ^-|6|, and are of great importance in a number of scientific disciplines. 
For example, lipid bilayers are an integral component of cells, where they form the outer 
membrane and also play a role in transport processes jsl- The bilayer vesicles that form 
from block copolymers, on the other hand, are longer-lived and less permeable than their 

n Q 

lipid counterparts [5[ making them promising candidates as vehicles for drug delivery [7[. 

The encapsulation of hydrophobic molecules between the two leaflets of the bilayer has 
been a recurrent issue, and has been discussed in a variety of contexts, including drug 
carrier design, vesicle formation from inverse phase methods, and lipid distribution in cell 
membranes. In particular, it is often important to know whether the hydrophobic molecules 
form droplets or are spread more evenly throughout the bilayer center. For example, lipid 
cell membranes are often found to contain a second species of lipid js-lO], in surprisingly 



high concentrations The question then arises of whether this is due to the presence 

of small lipid domains within the bilayer, the existence of which has also been suggested to 
explain the formation of the lipid droplets seen in the center of cells [121, 13] • Observations 



of such bilayer domains have indeed been reported in the biophysics literatur e |l4j . [iSj, and 
these structures have recently been seen in molecular dynamics simulations However, 
their probable small size and the existence of other lipid domains in the cell complicates the 
interpretation of experiments, and the issue is not fully resolved jl^. 

The formation of oil droplets in amphiphilic bilayers is also a problem of current interest 
in microfiuidics, and has been observed in recent experiments on block copolymer systems 
17|. The aim of this research was to produce aqueous solutions of monodisperse 



vesicles from block copolymers in water-oil-water double emulsions by evaporating the oil. 
Here, the presence of oil droplets in the bilayers is undesirable, as it leads to unevenness in 
the vesicle wall. 



In addition, the encapsulation of hydrophobic mo 
delivery of drugs using block copolymer vesicles 120 



ecules in bilayers is of importance in the 



23|. Although it might at first appear 



more natural to encapsulate a hydrophobic substance in the core of a spherical micelle, 
copolymer vesicles can offer certain advantages over these smaller structures. In particular. 
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they can encapsulate both hydrophobic and hydrophihc compounds 24. |25 [. Furthermore. 



26|. 



faster release of the hydrophobic compound can be obtained from vesicles 

In this paper, we investigate how much information about hydrophobic droplets in bilayers 
can be obtained from a simple mean- field model of oil and amphiphile in solution. First, 
we will study the shape of the droplet and to what extent this varies with its size. We will 
then relate our results to capillarity arguments. Next, we will investigate the effect of the 
strength xbo of the repulsion between the hydrophobic liquid and the amphiphile tails on the 
droplet shape. This question is of relevance to several of the situations described above, and 
our results will give some guidance as to how robust the phenomenon of droplet formation 
is expected to be in experiments. In addition, understanding the role of the interaction 

n 

strength might allow an oil to be chosen to encourage or discourage [17J the formation of 
well-defined droplets. Finally, with similar objectives in mind, we will study the effect of 
the length of the hydrophobic molecules on the droplet shape and stability. 

The paper is organized as follows. In the following section, we introduce the theoretical 
technique to be used, self-consistent field theory. We then present and discuss our theoretical 
results, and give our conclusions in the final section. 



II. SELF-CONSISTENT FIELD THEORY 



Self-consistent field theory (SCFT) [27| has been used with success to investigate the 



equilibrium structures formed in melts and blends of polymers 28l430| . and may also be used 



to study metastable structures. 
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32l | and amphiphiles in solution |33 



34 



to a wide rarige of amphiphilic molecules, including simple homopolymers 
copolymers jssi, Is^] and any given mixture of these 



It can be applied 



35| , more complex 



SCFT requires less computational 



power than simulation methods such as Monte Carlo, yet often provides comparably accurate 
predictions of the form of individual structures 33|, |39|, |40|. Furthermore, as a coarse-grained 
model, with a simple description of the polymer molecules, it will allow us to capture the 
basic phenomenology of the system clearly. 



We now give a brief introduction to SCFT, and refer the reader to reviews 
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41 



42| for 



a fuller presentation. A complete description of our calculations for amphiphiles in solution 



is given in a recent publication 



431], and we present details only when our current system 



differs from that described there. SCFT models individual molecules as random walks in 
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space, and so neglects fine details of their structure and packing [42||. An ensemble of 
many such molecules is considered. The interactions between the molecules are modeled by 
assuming that the blend is incompressible and by introducing contact potentials between the 



molecules 



30l |. The stren gths of the potentials between the various species are specified by 



the Flory parameters Xij The computational difficulty of the problem is then sharply 
reduced by making a mean- field approximation [30]; that is, by neglectinp; fluctuations. 
This approximation is quantitatively accurate when the molecules are long 
addition, SCFT can provide considerab 
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4l|. In 



molecules, such as lipid bilayers 
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e qualitative insight into systems containing smaller 
43| and aqueous solutions of copolymer |34l. |45|. 
We now discuss the application of SCFT to our system of amphiphile and oil in a solvent, 
which we model by a mixture of block copolymer with two incompatible homopolymers that 



represent the oil and the solvent respectively. Although such a mixture of polymers may 



appear quite simple, models of this 



of hpid and copolymer systems [32 



evel of complexity have been used to study a wide range 



46| , and can capture broad phenomenology more clearly 



than more complicated theories. We take the copolymer to have a mean-squared end-to-end 
distance of a'^N, where a is the monomer length and is the degree of polymerization 30 |. 
One half of the monomers in this polymer are hydrophilic (type A) and the other half are 
hydrophobic (type B), so that the degrees of polymerization Aa and A^b for the A and B 
blocks are equal. We choose the same value of a'^N for the A homopolymer solvent as for 
the copolymers. Together with the values of and A^b, this ensures that the amphiphile 
preferentially forms flat bilayer structures [43] for the interaction strength we will consider 
here. The degree of polymerization A"o = aN of the oil will be varied between N/ 4 and 2A^. 

In this paper, we keep the amounts of copolymer and homopolymer in the simulation box 
flxed; that is, we work in the canonical ensemble. This will make it easier for us to access 
more complex structures such as droplets. Such structures are more difficult to stabilize in 
ensembles where the system is able to relax by varying the amount of the various species, 
and can require constraints to be imposed on the density profile j32 |. 

For concreteness and to introduce the appropriate notation, we note that the SCFT 
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approximation to the free energy of our system has the form 
FN _ Fi,N 

- {l/V) j dr [xiVAB(0A(r) + 0s(r) - a " 0s)(0B(r) - 0b) 

+ X^Ao(0A(r) + 0s (r) - 0A - 0s)(</'o(r) - 0o) 
+ xiVBo(0B(r) - 0B)(0o(r) - 0o)] 

- (0A + 0b) HQab/V) - 0s ln(gs/V^) - (0o/«) ln(go/V^) (1) 

where the 0j are the mean volume fractions of the various components. The 0j(r) are the 
local volume fractions, with i = A for the hydrophilic blocks, z = B for the hydrophobic 
blocks, i = O for the oil, and i = S for the solvent. The first Flory parameter, xab, is 
set to 50/A^, so that sharp, well-defined bilayers form. The other Flory parameters will 
be varied to study the effect of the nature of the oil on the droplet shape. V is the total 
volume of the system, 1/po is the volume of a monomer, and Fh is the SCFT free energy 
of a homogeneous system containing the same components. The details of the individual 
polymers enter through the single-chain partition functions Qi. These are calculated js^ 
from integrals over the propagators q and q\ which are also used to compute the polymer 



density profiles 
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41| . Reflecting the fact that the molecules are modeled as random 



Des 



walks, the propagators satisfy modified diffusion equations with a field term that descri 
the polymer interactions. These equations are solved using a finite difference method 
with step size of 0.04 aN^/"^. We assume that the droplet is cylindrically symmetric and forms 
at the center of the system, and hence consider an effectively two-dimensional problem in a 
cylindrical calculation box. Reflecting boundary conditions are imposed at the edges of the 
system. 

The derivation of the mean-field free energy F also generates a set of simultaneous equa- 
tions linking the values of the fields and densities. In order to calculate the SCFT density 
profiles for a given set of polymer concentrations, we begin by making an initial guess for 
the fields Wi(r) and solve the diffusion equations to calculate the propagators and then the 
densities corresponding to these fields. The new 0i(r) are then substituted into the simul- 
taneous equations to calculate new values for the Wi j48|. The procedure is repeated until 
convergence is achieved. We have checked that the algorithm converges to the same solution 
from different initial states and with different iteration speeds. 



To form the structure we wish to study, a bilayer in the z = plane encapsulating a 
droplet at its center, a suitable initial guess for the Wi must be made. It is important to 
note that it is not necessary to include any detailed information about the shape of the 
droplet or bilayer in this ansatz. Although, for the sake of speed, we will often use the 
final self-consistent fields corresponding to one bilayer-droplet system as initial guesses for 
a subsequent calculation, the initial form of the fields can be very simple. Specifically, to 
encourage the formation of the bilayer, it is sufficient to start the SCFT iteration with a 
simple square well for the hydrophobic block field wb- If zq is the approximate width of 
the bilayer, we set wb to a low value for —z^jl < z < Zq/2, and a higher value elsewhere. 
The initial value of the hydrophilic block field wx can simply be set to zero, as the A and B 
blocks are connected and the above ansatz for wb is enough for an AB bilayer to form. The 
iteration for the field corresponding to the oil, wq, can similarly be initiated with a square 
well potential. The potential is set to a low value in a cylindrical region at the center of the 
system {—z'q < z < z'q, r < tq), with z'q < zq and vq much smaller than the radius of the 
simulation box. 

We set the tension of a bilayer with no oil to zero, as this corresponds most closely 
to the experimental situation of a vesicle in solution. To find the zero-tension bilayer, we 
proceed as follows 46] . First, we calculate the free-energy density of a (one-dimensional) box 
containing an infinite planar aggregate in solvent. The volume of the simulation box is then 
varied in the direction perpendicular to the bilayer surface, keeping the volume fraction of 
copolymer fixed (at 10%), until the box size with the minimum free-energy density is found. 
This scheme was introduced to mimic the behavior of a system of many aggregates 46, 49l |. 
which minimizes its free energy by varying the number of aggregates and hence the volume 
('box size') occupied by each. Equivalently, this procedure allows us to prepare a bilayer 
under zero tension. Decreasing the box size at constant copolymer volume fraction reduces 
the amount of amphiphile in the system and so thins the bilayer, which corresponds to 
stretching it parallel to its surface. Conversely, increasing the box size thickens the bilayer, 
which is physically equivalent to compressing it. The bilayer found for the box size where 
the free energy is at a minimum is one that is neither too stretched nor too compressed and 
has no contributions to the free energy from polymer chains that are forced into unfavorable 
configurations. 

This calculation is then used to fix the size of the cylindrical box in the 2;-direction 
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to 17.6 aA^^/^, so that —Z < z < Z, where Z = 8.8. The radius of the box is set to 
R = 16 aN^^"^^ to allow droplets of a wide range of sizes to be studied. 

III. RESULTS AND DISCUSSION 

In this section, we study the shape and size of the droplets in detail for a single set of 
system parameters. Then we relate the droplet shape to the bilayer and monolayer tensions 
and amphiphile concentrations. Finally, we investigate how the droplet shape and stability 
depends on the nature and size of the oil molecules. 

A. Droplet shape. 

To begin, we calculate the density profiles for bilayer-encapsulated droplets of various 
sizes. We focus on a system with oil molecules that are half the size of the amphiphiles, 
so that Nq = N/2. The interaction strength between the hydrophobic B-block of the 
amphiphile and the oil is set to xbo = The effect of varying this quantity will be 

investigated later. Given that the strength of the repulsion between the A and B blocks of 
the amphiphile has already been set to xab = 50/A^, we no longer have complete freedom 
in our choice of the final Flory parameter, xao- If we assume that Xij is related to the 
polarizabilities ajj of the two polymer species by Xij = ~Q^j)^) where k is a constant of 



proportionality [50] , we find that xao is given in terms of the other two interaction strengths 



and is therefore set to xao = 23 A/ N, where we choose the negative root to give a moderate 
incompatibility between the oil and the solvent. 

Starting with the tensionless bilayer described above, we compute the density profiles 
of bilayer-encapsulated droplets for a range of oil volume fractions between (pQ = 0.01 
and 0o — O-l- The first of these values corresponds to the smallest droplet that could 
be stabilized in our calculations. In Figure [H we show the density profiles for (a) the 
amphiphile and (b) the oil for an intermediate-sized droplet with (pQ = 0.04. Figure [IK 
clearly shows the splitting of the amphiphile bilayer into two thin monolayers to incorporate 
the droplet. In Figure [T]d, we plot the density profile of this lens-shaped droplet, and see 




by 




(2) 
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also that a significant amount of oil remains between the two leaflets of the bilayer in the 
region surrounding the drop. This feature appears as horizontal gray lines on either side of 
the droplet in the density plot of Figure [T]d, and is a result of the relatively weak repulsion 
Xbo between the hydrophobic block and the oil. 

To help visualize the droplets, we show a ray-traced plot of the droplet surface (defined 
as the locus of points where 0o(r) = 0.5) in Figure |2l The overall lens shape of the droplet 
is clearly visible, as is the slight rim where the edge of the droplet meets the oil remaining 
in the bilayer. 



B. Surface concentrations and tensions. 



In order to gain more detailed insight into the droplet shapes, we plot cuts through the 
droplet surface for a range of droplet sizes from the smallest to the largest (see Figure |3]). 
If we assume that no long-range forces act in the system, that the bending rigidity of the 
membrane can be neglected, and that the pressure inside the droplet is constant, both the 



upper and lower 



Laplace law 



52 



lalves of the droplet will have a spherical cap shape 5l|] in order to obey the 



53| . The expected shape for these cuts shown in Figure |3] is then a section 



of a circle, with the contact angle 9 (see inset) determined by the mechanical equilibrium 
of the surface tensions along the contact line [Sjj]. We indeed find that this shape gives a 
very good fit to the data for all droplet shapes (Figure E]), with a slight deviation in the rim 
region shown in Figure |2l where the droplet spreads out slightly between the two amphiphile 
leaflets instead of forming a perfect cusp. In addition, the contact angle calculated from the 
fits is the same for all droplet sizes to the accuracy of the calculations, and is given by 
e ^ 51°. 

To understand how the presence of the bilayer leads to the formation of these lens-shaped 
droplets, and to illustrate some other features of the density profiles plotted in Figure [H we 
now plot a series of cuts through the density profiles of the various species. In Figure H^, 
we show a cut in the z-direction (perpendicular to the bilayer) through the density profiles 
of all species at the edge of the system containing the smallest droplet studied. This plot 
shows the state of the bilayer as far away from the droplet as possible. In Figure Hb, we 
show the corresponding plot for the bilayer in the system with the largest droplet. First, 
we note that the two plots are very similar and that the bilayer is not strongly distorted 
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by the presence of the droplet. We will return to this point in a more quantitative fashion 
later on. Second, we see that a small but significant amount of oil remains in the center of 
the bilayer. Furthermore, the density profile of the oil shows a clear peak at the interface 
between the hydrophilic A and hydrophobic B blocks. This is because both xao and xbo 
are smaller than xab, the repulsion between the two blocks of the amphiphile. A thin film 
of oil therefore forms in this region to protect these two strongly incompatible species from 
each other. 

A clear contrast is seen between the density profiles of the bilayer a long way from the 
droplet, and those of the monolayer that covers the droplet. These latter profiles, calculated 
at the center of the system (r = 0), are plotted in Figure Ht (for the smallest droplet) 
and d (for the largest droplet). We note that, in both these cases, the maximum values 
of the density profiles of the A and B blocks of the amphiphile are lower than in Figure 

and b. This is because the monolayer has to be stretched and thinned to cover the 
droplet, lowering the surface amphiphile concentration. In addition, the difference between 
the profiles for the smallest and largest droplets is much more marked than in the case of the 
bilayer. In particular, the amphiphile concentration in the monolayer covering the largest 
droplet (Figure HJi) is noticeably lower than in that covering the smallest droplet (Figure HI;), 
showing that the monolayer must be further stretched to encapsulate more oil. Furthermore, 
the monolayer in Figure Hll is more symmetric with respect to its inner and outer leaflets than 
that in Figure Hb, and the peak values of the A- and B-block concentrations are much closer. 
The reason for this is that the surface of the larger droplet is flatter, and the monolayer that 
encloses it is quite close to that which would form at a planar oil-solvent interface. 

We now present a more quantitative discussion of the amphiphile density profiles. To 
begin, we integrate the bilayer and monolayer density profiles (Figure H]) in the z-direction, 
including both hydrophilic and hydrophobic blocks. Specifically, we calculate 



Fi is then the surface density of (half) the bilayer at the edge of the system {r = R), while 
r2 is the monolayer surface density at the center of the system (r = 0). The amount of 
amphiphile remaining in the bulk is low, so that its contribution to the surface densities is 
very small. Fi and r2 are calculated for all droplet sizes studied and are plotted against 




(3) 
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each other in Figure |5^. As would be expected from the profiles shown in Figure HJ the 
monolayer surface density varies over a wider range than the corresponding quantity for the 
bilayer, as the amphiphilic molecules are spread out more and more thinly in the monolayer 
as the droplet size increases. Furthermore, Fi and F2 are linearly related for a wide range of 
droplet sizes, with deviations from linearity only being seen for the higher surface densities 
corresponding to very small droplets. To understand this behavior, we relate Fi and F2 to 
the corresponding surface energies. First, we note that the amphiphile in our system acts as 
a surfactant, separating the solvent from the oil in the droplet and bilayer. Then, we assume 
that adding amphiphile linearly reduces the surface tension from its value in the absence of 
amphiphile, 70, so that 7i = 70 — 5Fj. Here, 5 is a constant of proportionality and i = 1 
for the tensions and densities at the system edge and 2 for those at the droplet surface. 
Balancing these two tensions at the rim where the droplet meets the bilayer, as shown in 
the inset to Figure [3], we find also that 71 = 72 cos 0. Combining this with our expressions 
for 7i and 72, we have 



From the slope of the straight line in FigurelSH, we find that cos 6 ^ 0.627, so that 6 ^ 51°, 
in excellent agreement with our independent measurement of 6 from the cross-sections in 
Figure [3l This shows the validity of the force balance argument, and also confirms our use 
of the same proportionality constant 6 in our expressions for 71 and 72. To check our linear 
formula for 71, we have also calculated, using Equation [H the free energy density of a flat 
oil-containing bilayer with the profile shown in Figure H^. Similar calculations are performed 
for all values of 0o- W^e then plot the quantity / = FN/k-QTpoV — F^N/k-BTpoV (the free 
energy density measured with respect to that of the homogeneous solution with the same 
composition) as a function of the surface concentration Fi. As can be seen from Figure [5b, 
/ decreases linearly with Fi for all but the very smallest droplets, confirming our simple 
model for the surface energy 71. We note that the free energy density as calculated from 
Equation [1] includes a contribution from the solvent region as well as from the bilayer itself. 
However, this is likely to have a relatively small effect on the variation of /, as the bulk 
amphiphile concentration changes very little with droplet size. 

Both plots show some deviation from linearity for the smallest two or three droplet sizes 
considered. For the plot of the two surface concentrations in Figure [5^, the deviation comes 
from the increasing relative importance of the rim around the edge of the droplet (Figure 




(4) 
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[2]), which means that a simple force balance argument based on a well-defined contact angle 
is less valid. The slight breakdown of linearity in Figure [5]d may be due to the fact that the 
droplet is nearing its lower size limit and the free energy density / is dropping more rapidly 
as the bilayer relaxes towards the flat state. 



C. Nature and size of the oil molecules. 



Having established a basic picture of droplet formation in our system, we now turn our 
attention to the question of how the nature and size of the oil molecules affects the shape 
and stability of the droplets. First, we investigate the effect of changing the Flory parameter 
Xbo determining the strength of the interaction between the oil and the hydrophobic blocks 
of the amphiphile. This corresponds to changing the chemical nature of the oil while keeping 
the length of the oil molecules constant. We keep the same type of amphiphiles as used in 
the preceding section, so that and xab are unchanged. However, as noted earlier, the 
three Flory parameters cannot be varied completely independently [50|], and xao must be 
recalculated according to Equation [2] for each value of Xbo- The volume fraction of oil 
is set to (pQ = 0.04. In Figure [6], we plot the outlines of the droplets formed as Xbo^ 
is decreased from 5 (the value used in our earlier calculations) to 1 in steps of 1. The 
droplet with Xbo^ = 5 is that with the most rounded shape and the greatest thickness 
in the z-direction. As Xbo^ is lowered to 4, the droplet spreads outward slightly into the 
bilayer and becomes thinner, as the oil becomes more compatible with the hydrophobic 
blocks of the amphiphile. The spreading effect here is rather small, suggesting that, above a 
certain Xbo^; the droplet shape is relatively insensitive to the nature of the oil and retains 
its characteristic lens form. As Xbo^ is reduced further, the droplet continues to spread. 
However, the amount by which the droplet thickness falls due to each reduction of Xbo^ by 
1 gradually increases, until we obtain an almost flat structure at Xbo^ = 1- It is interesting 
to note that, even for this very weak repulsion between the oil and the hydrophobic blocks, 
the droplet remains at least metastable. 

The results shown in Figure [6] differ somewhat from the classical problem of the spreading 
of a single droplet 5l|], as our droplet is in equilibrium with a fllm of oil in the bilayer, which 
grows in thickness, taking material from the droplet, as the oil becomes less incompatible 
with the hydrophobic B-blocks. This is particularly clear for the lowest value of Xbo-^ 
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considered. Here, the oil concentration in the bilayer is so high that our definition of the 
droplet surface as the locus of points at which 0o(r) = 0.5 now includes the oil film as an 
extension of the droplet. 

Finally, we study the effect of the size of the oil molecules on the droplet shape. Returning 
to our original set of Flory parameters, with Xbo^ = 5, we consider the following values of 
the oil polymerization index: A^^o = 2A^, N, 0.5N (the original value), and 0.25A^, and plot 
the droplet outlines in Figure [71 As above, (pQ = 0.04. For the largest three values of Nq, 
the droplet shape changes rather little. It simply shrinks slightly as Nq is lowered, as this 
change reduces the repulsion between the oil and the hydrophobic sections of the amphiphile 
so that more material leaks out of the droplet into the oil film. 

However, as Nq is lowered further, to 0.25A^, a sharp change occurs in the droplet shape. 
The thickness of the droplet in the z-direction is now significantly greater, while its radius is 
smaller. It is difficult to interpret this result in terms of the simple force balance arguments 
used earlier. This is because the rim feature, which was previously a small perturbation on 
droplets whose shape could be represented by two joined spherical caps [51], is now a much 
more significant part of the droplet, since the small oil molecules penetrate more effectively 
into the bilayer. We speculate that the elongated shape of this droplet may be a precursor 
to its eventual splitting into two smaller monolayer-wrapped droplets, such as those recently 
studied by Kusumaatmaja and Lipowsky in the context of membranes in contact with several 
fluids 5J]. In any case, it certainly seems that the simple single-droplet solution to SCFT 
becomes unstable around Nq = 0.25N. If we reduce Nq below this value, no solution to the 
SCFT equations can be found using our current methods, and our algorithm slows down 
considerably even for A^o = 0.25A^. 



IV. CONCLUSIONS 



Using a coarse-grained mean-fleld approach (self-consistent fleld theory) we have modeled 
several aspects of the structure of hydrophobic droplets encapsulated between the two leaflets 
of an amphilic bilayer. First, we have found that droplets of a range of sizes have the same 
simple lens shape that would be expected from simple capillarity arguments. We have 
explained both this shape and the amphiphile concentrations in different regions of the 
system by considering the balance of the surface tensions around the edge of the droplet. 
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Next, we studied the effect of the oil parameters on the droplet shape. We found that, 
although reducing the incompatibility xbo causes the droplet to flatten and spread outwards 
into the bilayer, it remains at least metastable even for very low Xbo- There appears to be 
no clear threshold value of xbo for droplet solutions to SCFT to exist. The droplets are also 
relatively insensitive to changes in the oil molecule length. In fact, provided these molecules 
are longer than the hydrophobic part of the amphiphile, their length has little effect on 
the droplet shape. These observations provide some evidence that droplet formation is 
a relatively robust effect and may therefore be a reasonable explanation for phenomena 



such as the inclusion of significant amounts of a second lipid species in lipid bilayers 
Furthermore, the formation of droplets even in our simple model indicates that this might 
be quite a general phenomenon and that hydrophobic domains such as those seen in the 



might be observed in a variety of 



molecular dynamics simulations of Khandelia et al. 
systems. 

The current work opens a number of interesting perspectives that could be discussed 
within the framework of self-consistent field theory. First, our free energy calculations 
could be extended, to find the parameter range where the oil will form a droplet rather 
than spreading. Second, the question of whether an optimum droplet size exists could be 
addressed, perhaps by using a range of system sizes or by considering the stability of a 
droplet with respect to two smaller droplets. We could also study how likely the droplet 
is to split off from the bilayer, for example by comparing the free energies of the bilayer- 
encapsulated droplet and a system of a droplet covered by a monolayer in coexistence with 



a bilayer. Finally, droplets formed from a second species of amphiphile coulc 
bring our calculations closer to the problem of lipid domains in bilayers 
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FIG. 1: Density plots of (a) amphiphile (including both hydrophilic and hydrophobic blocks) and 
(b) oil in a droplet-bilayer system with (j)Q = 0.04. Cylindrical polar coordinates are used, and 
dark regions indicate high volume fraction. 




FIG. 2: Ray-traced plot of the surface of the droplet shown in Figure [T]3. 




FIG. 3: Cuts through the droplet surface for a range of droplet sizes with <j)Q = 0.01, 0.02, 0.04, 
0.06, 0.08, and 0.1. Circles show the data points from our SCFT calculations. Although all points 
are used in the fits, only every fifth point is plotted for clarity. Solid lines show fits to the data 
using sections of a circle. The inset shows the contact angle 9 and the surface tensions 71 and 72, 
for later reference. 
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FIG. 4: Cuts through the density profiles perpendicular to the bilayer at the edge of the system 
and at the center of the droplet. The hydrophilic A-blocks are shown with thick dashed lines, 
the hydrophobic B-blocks with thick full lines, the oil with thin full lines, and the solvent with 
thin dashed lines, (a) Bilayer at the edge of a system containing the smallest droplet studied 
(^O = 0.01). (b) Bilayer at the edge of a system containing the largest droplet (^o = 0.1). (c) 
Monolayer covering a droplet at the center of the = 0.01 system, (d) Monolayer covering a 
droplet a the center of the ^q = Q.\ system. Note the change in z-axis between (c) and (d). 
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FIG. 5: (a) Surface density of the monolayer covering the droplet plotted against that of (half) the 
bilayer at the edge of the system, (b) Normalized free energy density of the bilayer plotted against 
its surface density. 
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FIG. 6: Droplet outline for a range of values of Xbo^ from 1 (flattest drop) to 5 (roundest drop). 

The oil volume fraction 0o is fixed to 0.04. 




FIG. 7: Droplet outline for a range of values of Nq from 2N (outermost flat drop), through N and 
0.5N, to 0.25iV (round drop). 
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